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ABSTRACT. Multi-frequency, high resolution, full sky measurements of the 
anisotropy in both temperature and polarisation of the cosmic microwave 
background radiation are the goals of the satellite missions MAP (NASA) 
and Planck (ESA). The ultimate data products of these missions — multiple 
microwave sky maps, each of which will have to comprise more than ~ 10 6 
pixels in order to render the angular resolution of the instruments — will 
present serious challenges to those involved in the analysis and scientific ex- 
ploitation of the results of both surveys. Some considerations of the relevant 
aspects of the mathematical structure of future CMB data sets are presented 
in this contribution. 



1 Introduction 

The extraordinary success of 
NASA's COBE satellite mission (see 
www.gsfc. nasa. gov / astro / cobe) , and 
in particular the COBE-DMR dis- 
covery of anisotropy in the temper- 
ature of the cosmic microwave back- 
ground (CMB) radiation (Smoot et 
al. 1992), has transformed the field 
of experimental and theoretical stud- 
ies of the CMB as an astronom- 
ical window on the early universe 
into one of the most vibrant areas 
of modern cosmology. The explosive 
growth of activity in the CMB com- 
munity revealed the promise of dra- 
matic improvement in our knowl- 
edge of the early universe that 
should materialize given the contin- 
uation of efforts to study the CMB 
from space. This, and the realization 
of a sufficient improvement in mi- 
crowave detector technology suitable 
for space experiments, resulted in 
both NASA and the European Space 
Agency approving proposals to con- 
duct the next generation satellite 
miss ions to study the CMB MAP 
(see |http:/ /map. gsfc.nasa.gov ) to be 
launched in 2000. and Plan ck (see 



fittn: / /astro, cstcc.esa.nl /SA- 



to 



be 



jeiieral/Projects/Planck ) 
launched in 2007. 

Both MAP and Planck were de- 
signed to surpass the COBE-DMR 
capabilities to measure anisotropy in 
the temperature and polarization of 
the CMB by orders of magnitude. 
Both missions will observe the CMB 



from the Earth-Sun L-2 point us- 
ing ~ 1.5m telescopes. MAP will 
use passively cooled HEMT detec- 
tors with five frequency bands be- 
tween 22 and 90 GHz, and will reach 
an angular resolution between 0.93 
and 0.21 deg, respectively, and is ex- 
pected to render a noise sensitiv- 
ity of ~ 20a*K per 0.3 x 0.3 de- 
gree pixel achieved by combining 
the three highest frequency chan- 
nels. Planck will comprise two instru- 
ments involving complementary de- 
tector technologies to yield unprece- 
dented, very wide frequency cover- 
age. The Low Frequency Instrument 
will use an actively cooled, tuned ra- 
dio receiver array operated at 20 K 
with four frequency bands between 
30 and 100 GHz, should achieve an 
angular resolution between 33 and 10 
arcmin, respectively, and should ren- 
der a noise sensitivity per resolution 
element between 4 and 12 /iK, re- 
spectively. The High Frequency In- 
strument will use bolometer arrays 
operated at 0.1 K with six frequency 
bands between 100 and 857 GHz, 
should achieve an angular resolution 
between 10 and 5 arcmin, respec- 
tively, and should render a noise sen- 
sitivity per resolution element be- 
tween 4.6 and 12 fj,K in the CMB 
dominated channels (100-217 GHz). 

It is widely recognized that if both 
missions succeed in achieving these 
spectacular specifications during ob- 
servation of the CMB sky our knowl- 
edge of the universe will be fur- 
thered dramatically (see the con- 
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tributions by C. Lawrence and N. 
Sugiyama, and references therein, in 
these proceedings). Indeed, it is ex- 
pected that the results of MAP and 
Planck should allow us to answer de- 
cisively both fundamental questions 
about the global properties of the 
universe — the average density of 
its matter content, the expansion 
rate, global curvature, the existence 
of a cosmological constant and/or a 
cosmological background of gravity 
waves, to name a few — and to un- 
ravel the properties of a number of 
astronomical objects — for exam- 
ple the detailed picture of emission 
from our Galaxy, and infrared emis- 
sion from distant galaxies and galaxy 
clusters. 

At the same time it should be re- 
alized that the data products of the 
future CMB missions will be by no 
means trivial to work with. The enor- 
mity of the expected scientific return 
will come at the price of a concerted 
effort to improve all currently avail- 
able methods of analysis of very large 
CMB data sets, which, in order to 
answer properly some of the grand 
questions, will need to be studied 
globally. 

As we have learned working with 
the COBE mission products, the 
digitized sky map is an essential 
intermediate stage in information 
processing between the entry point 
of data acquisition by the instru- 
ments — very large time ordered 
data streams, and the final stage 
of astrophysical analysis — typically 
producing a 'few' numerical values 
of physical parameters of interest. 
COBE-DMR sky maps (three fre- 
quency bands, two channels each, 
6144 pixels per map) were consid- 
ered large at the time of their re- 
lease. So, looking forward to the gi- 
ant leap to be taken by MAP and 
Planck in mapping the entire CMB 
sky, what do we mean when saying 
'very large future CMB data sets'? 
Very simply, if we focus on the case 
of whole sky CMB measurements at 
the angular resolution of ~ 10 ar- 
cmin (FWHM), and presume that 
a few pixels per resolution element 
should be used to discretise the sig- 
nal in a non-damaging way (i.e. in 
such a way that discretisation effects 
are sufficiently sub-dominant with 
respect to the effects of instrument 
response), we require that the entire 



map should comprise at least N p i x ~ 
a few X 1.5 10 s pixels. More pixels 
than that will be needed to repre- 
sent the Planck-HFI higher resolu- 
tion channels. This estimate, N p i x , 
should then be multiplied by the 
number of frequency bands (or, in- 
deed, by the number of individual 
observing channels — 74 in the case 
of Planck — for the analysis work 
to be done before the final coadded 
maps are made for each frequency 
band) to render an approximate ex- 
pected size of the already compressed 
form of data, which would be the 
input to the astrophysical analysis 
pipeline. Clearly, it is easy to end 
up with an estimate of many GBy. 
Hence it appears that some atten- 
tion should be given to devising such 
high resolution CMB map structures 
which would maximally facilitate fu- 
ture analyses of large size data sets. 

The main part of this contribu- 
tion is devoted to a presentation of 
the properties of one such proposed 
approach to high resolution full sky 
map making — the Hierarchical 
Equal Area and isoLatitutde Pix- 
elizat ion of the sphere (HEALPIX . 



http:/ /www.tac.dk/^hcalpix) 



The remaining part of this paper is 
a digression on the issue of the (non- 
)gaussianity of the CMB sky. 



2 Discretisation of Functions 
on a Sphere with HEALPIX 
for High Resolution Appli- 
cations 

Operations which are, or will be, 
routinely executed in the analy- 
sis of CMB maps include convo- 
lutions with local and global ker- 
nels, Fourier analysis with spheri- 
cal harmonics and power spectrum 
estimation, wavelet decomposition, 
nearest-neighbour searches, topolog- 
ical analysis (minimum/maximum 
search, Minkowski functionals, etc.), 
and more. Some of these operations 
become prohibitively slow if the sam- 
pling of functions on a sphere, and 
the related structure of the discrete 
data set, is not designed carefully. 

Typically, a whole sky map ren- 
dered by a CMB experiment contains 
the sky signals, which should be (by 
design) strongly band-width limited 
(in the sense of spatial Fourier de- 
compostion) by the instrument's re- 
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Figure 1. Orthographic view of HEALPIX division of a sphere. Overplot 
of equator and meridians illustrates octahedral symmetry of the HEALPIX 
construction. Light-gray shading shows one of eight (four north, and four 
south) identical polar base-resolution pixels. Dark-gray shading shows one 
of four identical equatorial base-resolution pixels. Moving clockwise from the 
upper left panel the grid is hierarchically subdivided with the grid resolution 
parameter equal to N ail { e = 1, 2, 4, 8, and the total number of pixels equal 
to N pix = 12 x N^ ide = 12, 48, 192, 768. All pixel centers are located on 
Nring = 4 X N ai ^ e — 1 rings of constant latitude. Within each panel the areas 
of all pixels are identical. 



sponse function, and projection of 
the instrumental noise, which, at 
least near the discretisation scale, 
should be random, with a band- 
width significantly exceeding that of 
all the signals. 

With all these considerations in 
mind one may propose the following 
list of desiderata for the mathemat- 
ical structure of discrete whole sky 
maps: 

1. Hierarchical structure of 
the data base. This is recog- 
nized as essential for very large 
data bases, and was indeed postu- 
lated already in construction of the 
Quad rilatcralizcd Spherical Cube 
(see fittp: / / www.Ksfc.iiasa.gov/ast- 
ro / cobc/sky map_info. htird] ) , which 
was used for all the COBE data. 
A simple argument in favour of this 



states that the data elements which 
are nearby in a multi-dimensional 
configuration space (here, on the sur- 
face of a sphere), are also nearby 
in the tree structure of the data 
base. This property facilitates vari- 
ous topological methods of analysis, 
and allows for easy construction of 
wavelet transforms on triangular and 
quadrilateral grids. 

2. Equal areas of discrete el- 
ements of partition. This is ad- 
vantageous because white noise at 
the sampling frequency of the in- 
strument gets integrated exactly into 
white noise in the pixel space, and 
sky signals are sampled without re- 
gional dependence (although still 
care must be taken to choose a pixel 
size sufficiently small compared to 
the instrumental resolution to avoid 
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excessive, and pixel shape depen- 
dent, signal smoothing). 

3. Iso-Latitude distribution of 
discrete area elements on a 
sphere. This property is essential 
for computational speed in all opera- 
tions involving evaluations of spher- 
ical harmonics. Since the associated 
Lcgcndrc polynomials are evaluated 
via slow recursions, any sampling 
grid deviations from an iso-latitudc 
distribution result in a prohibitive 
loss of computational performance 
with the growing number of sampling 
points. 

Various known sampling distribu- 
tions on a sphere fail to meet simul- 
taneously all of these requirements: 

i) Quadrilatcralized Spherical 
Cube obeys points 1 and (approxi- 
mately) 2, but fails on point 3, and 
cannot be used for efficient Fourier 
analysis at high resolution. 

ii) Equidistant Cylindrical Projec- 
tion, a very common computational 
tool in geophysics, and climate mod- 
cling, and recently implemented for 
CMB work (Muciaccia, Natoli, Vit- 
torio, 1998), satisfies points 1 and 3, 
but by construction fails with point 
2. This is a nuisance from the point 
of view of application to full sky sur- 
vey data due to wasteful oversam- 
pling near the poles of the map. 

iii) Hexagonal sampling grids with 
icosahcdral symmetry perform su- 
perbly in those applications where 
near uniformity of sampling on a 
sphere is essential (Saff, Kuijlaars, 
1997), and can be devised to meet re- 
quirement 2 (see e.g. Tegmark 1996). 
However, by construction they fail 
both points 1 and 3. 

iv) Igloo-type constructions are 
devised to satisfy point 3 (E. Wright, 
1997, private communication; Crit- 
tenden & Turok, 1998). Point 2 
can be satisfied to reasonable ac- 
curacy if quite a large number of 
base-resolution pixels is used, which 
precludes the efficient construction 
of simple wavelet transforms. Con- 
versely, a tree-structure seeded with 
a small number of base-resolution 
pixels forces significant variations in 
both area and shape of the pixels. 

All three requirements formulated 
above are satis- 

fied by construction with the Hierar- 
chical Equal Area, iso-Latitude Pix- 
elisation (HEALPIX) of the sphere 
(Gorski 1998, Hivon, Gorski, 1998), 



which is shown in Figure 1. 

HEALPIX is a genuinely curvi- 
linear partition of the sphere into 
exactly equal area quadrilaterals of 
varying shape. The base-resolution 
comprises twelve pixels in three rings 
around the poles and equator. 

The resolution of the grid is ex- 
pressed by parameter N^,, which 
defines the number of divisions along 
the side of a base-resolution pixel 
that is needed to reach a desired 
high-resolution partition. 

All pixel centers are placed on 
rings of constant latitude, and are 
equidistant in azimuth (on each 
ring). All iso-latitude rings located 
between the upper and lower corners 
of the equatorial base-resolution pix- 
els, or in the equatorial zone, are di- 
vided into the same number of pix- 
els: N eq = 4 X -/Vside- The remaining 
rings are located within the polar cap 
regions and contain a varying num- 
ber of pixels, increasing from ring to 
ring, with increasing distance from 
the poles, by one pixel within each 
quadrant. 

Pixel boundaries are non-geodesic 
and take the very simple form: 
cos 8 = a + b X (j> in the equato- 
rial zone, and cos = a + b/(f> 2 in 
the polar caps. This allows one to 
explicitly check by simple analyti- 
cal integration the exact area equal- 
ity among pixels, and to compute ef- 
ficiently more complex objects, e.g. 
the Fourier transforms of individual 
pixels. 

Specific geometrical properties al- 
low HEALPIX to support two dif- 
ferent numbering schemes for the 
pixels, as illustrated in the Figure 
2. First, one can simply count the 
pixels moving down from the north 
to the south pole along each iso- 
latitude ring. It is in this scheme 
that Fourier transforms with spher- 
ical harmonics are easy to imple- 
ment. Second, one can replicate the 
tree structure of pixel numbering 
used e.g. with the Quadrilatcralized 
Spherical Cube. This can easily be 
implemented since, due to the simple 
description of pixel boundaries, the 
analytical mapping of the HEALPIX 
base-resolution elements (curvilin- 
ear quadrilaterals) into a [0,1] X [0,1] 
square exists. This tree structure, 
a.k.a. nested scheme, allows one to 
implement efficiently all applications 
involving nearest-neighbour searches 
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Figure 2. Cylindrical projection of the HEALPIX division of a sphere and 
two natural pixel numbering schemes (ring and nested) allowed by HEALPIX. 
Both numbering schemes map the two dimensional distribution of discrete 
area elements on a sphere into the one dimensional, integer pixel number 
array, which is essential for computations involving data sets with very large 
total pixel numbers. $From top to bottom: Panel one (resolution parame- 
ter N ail i e = 2) and panel two (N ail i e =4) show the ring scheme for pixel 
numbering, with the pixel number winding down from north to south pole 
through the consecutive isolatitude rings. Panel three (resolution parameter 
N s ide = 2 J and panel four (Ngi^ = A) show the nested scheme for pixel 
numbering within which the pixel number grows with consecutive hierarchical 
subdivisions on a tree structure seeded by the twelve base-resolution pixels. 



(see Wandelt, Hivon, and Gorski 
1998), and also allows for an imme- 
diate construction of the fast Haar 
wavelet transform on HEALPIX. 



We have developed a HEALPIX 
software package which contains a 
suite of programs to simulate and 
analyse full sky CMB maps, display 
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the results, manipulate the trans- 
formations from physical space to 
pixel number (and the reverse), in 
both numbering schemes, and con- 
duct nearest-neighbour and max- 
ima/minima searches on the maps. 
The p ackage is available to the pub- 
lic at littp://www. tac.dk/~healpix. 
HEALPIX is presently the format 
chosen by the MAP collaboration to 
be used for the production of sky 
from the mission data (see 





tittD: / /map. gsfc.nasa.gov /html/tech- 


nical_info.html). HEALPIX software 



is widely used for simulation work 
within both LFI and HFI consortia 
of Planck collaboration. 



in the DMR data analysis — and 
again, quite unlike the original effect 
in the multipole index dependence. 

Hence, our (FMG) original con- 
clusions remain unaffected: unless an 
as yet completely unrecognized fore- 
ground emission results in the de- 
tected non-Gaussianity, the remain- 
ing candidate for the source of this 
effect is the CMB anisotropy itself. 
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3 Digression: Non- 
Gaussianity in the DMR 
data? 

The remainder of this contribution 
is a digression, which was stimulated 
by questions during the presentation. 

In a recent paper (Ferreira, 
Magueijo, and Gorski, 1998) the au- 
thors argued that an estimate of 
the reduced bi-spectrum of the CMB 
anisotropy in the COBB-BMR 53 
and 90 GHz 4yr data is significantly 
deviant from the distribution sim- 
ulated from a Gaussian ensemble. 
Their work addressed the possibility 
of the effect being due to foreground 
emission from our Galaxy, but failed 
to explain it as originating from the 
known distributions of galactic emis- 
sion at low and high frequencies (as 
traced by radio synchrotron and dust 
maps respectively). 

An update on the status of our 
(FMG) ongoing tests and searches 
for alternatives tends to the con- 
clusion that the CMB itself may 
have revealed its non-Gaussian na- 
ture. Evidence in support of this is as 
follows: 1) the reduced bi-spectrum 
estimated on the DMR difference 
maps (noise-only combinations of the 
DMR sky maps) shows no excess at 
the conspicuous multipole index of 
I = 16 where the dominant contribu- 
tion to the original effect was found; 
2) likewise, the reduced bi-spectrum 
estimated on a suite of systematic ef- 
fect model templates that were de- 
rived during the CO-RE-DMR mis- 
sion, turned out to be insignificant 
in amplitude — itself a testimony to 
the superb systematic effect control 
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